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The paper is motivated by observation of a kind of branes in the vacuum state 
of the lattice SU{2) gluodynamics. The branes represent two-dimensional vortices 
whose total area scales in physical units while the non-Abelian action diverges in 
the ultraviolet. We consider the question whether effects of the branes can be 
accommodated into the continuum theory. We demonstrate that at least in case 
of the gluon condensate (plaquette action) and of the heavy quark potential the 
contribution of the branes corresponds to the ultraviolet renormalon. Thus, the 
vortices might represent a non-perturbative match of the ultraviolet renormalon. 
Such an identification constrains, in turn, properties of the branes. 



1. Introduction 

Recently it has been discovered that monopoles and central vortices- 
non-perturbative fluctuations commonly considered responsible for the 
confinement- have non-trivial structure in the ultraviolet. Namely both 
the monopoles, see ^ and references therein, and vortices ^ are associated 
with an excess of the non-Abelian action which is divergent in the ultravi- 
olet: 

L A 

{Sraon) ^ In 7 • — , {S^ort) ~ 0.54-^ , (1) 

where L is the length of the monopole trajectory, A is the area of the 
vortex, a is the lattice spacing representing the ultraviolet cut off. In case 
of monopoles, the overall constant, which we quote as In 7 is known actually 
with rather poor accuracy but this is not crucial for our purposes. 

Naively, one would expect that monopoles and vortices with action (1) 
propagate only very short distances, L ^ a, A r-^ a^. However, both 
monopoles and vortices form clusters which percolate through the whole 
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of the lattice volume V4. Defining the corresponding densities as: 

Lperc = ^PpercVi , Ayort = QpvortV4 , (2) 

one finds that both the monopole density (see, e.g., ^ and references therein) 
and the vortex density (for review and references see scale in physical 
units and are independent on the lattice spacing. According to the latest 
measurements: 

Pperc = 7.70(8) , A^ort « 4.0(2) /m-2 , (3) 

see ^ and ^, respectively. 

The data (1), (3) imply that the standard picture of vacuum fluctu- 
ations is to be adjusted to incorporate fine tuning, that is coexistence of 
the ultraviolet and infrared scales within the same fiuctuations A com- 
prehensive theory of the monopoles and vortices is not yet in sight mainly 
because the monopoles and vortices arc defined not in terms of the original 

Yang-Mills fields but rather in terms of projected fields, for review see, e.g., 
4 

While lattice data on the monopoles and vortices have been accumu- 
lating since long, the discovery of the ultraviolet divergences (1) makes the 
challenge to the theory much more direct. Indeed, it is commonly believed 
that in an asymptotically free theory all the ultraviolet divergences can be 
understood from first principles. And it is worth emphasizing at this point 
that the brane properties we are considering, that is (1), (3) are perfectly 
SU{2) invariant. 

We will argue that, indeed, starting from the continuum theory it is 

possible to derive strong constraints on the ultraviolet properties of non- 
perturbative fluctuations Moreover, the data can be confronted with the 
constraints without knowledge of the anatomy of the branes in terms of the 
original fields. 

The idea of the derivation is as follows. The continuum theory is well 
defined at short distances, while the branes are extended objects. Thus, 
we will project the eff'ect of the branes onto matrix elements calculable 
in the continuum theory. In particular, we will concentrate on the gluon 
condensate and heavy quark potential at short distances. 

Because of the asymptotic freedom the parton-model predictions are a 
reliable zero-order approximation for observables of this kind. If one tries 
to derive a complete answer, then the outcome is an infinite perturbative 



"we briefly mentioned the constraints in the talks 
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series plus power-like corrections. The power corrections arc the price for 
the factorial growth of the coefficients of perturbative expansions ^, for 
review see, e.g., 

We will argue that the brancs with the newly discovered properties do 
fit this scheme providing a match to the so called ultraviolet renormalon. It 
has never been foreseen, though, that a non-perturbative completion of the 
ultraviolet renormalon could be a brane with an action density divergent in 
the ultraviolet. Still, the new addition seems to fit well the phenomenology 
of the power corrections. In turn, identification of the branes with a non- 
perturbative counterpart of the ultraviolet renormalon puts constraints on 
properties of the brancs. 

Although our main interest is the effects of the branes, we begin in Sect. 
2 with a rather detailed discussion of the Q~^, or ultraviolet-renormalon 
related power corrections. In particular, we argue that there arc two com- 
plementary ways to account for these corrections. Either one confines one- 
self to low-order perturbation theory and adds the quadratic corrections 
ad hoc, or one is prepared to consider high orders in perturbation theory 
(sometimes, say, up to ten loops) and then the Q^^ corrections arc implic- 
itly contained in the perturbative sum. While we do not suggest any new 
explanation of a particular effect, our overall conclusions vary from what 
can be found in the literature. The point is that commonly it is assumed, 
explicitly or tacitly, that the two procedures outlined above are mutually 
inconsistent. 

In Sect. 3 we review phenomenological evidence on the corrections 
and argue that the data are consistent with the theoretical conclusions. 

In Sect. 4 we argue that the branes provide a non-perturbative match 
to the ultraviolet renormalon. In view of the discussion above this means, 
in turn, that the branes-when projected onto local matrix elements- are 
dual to high orders of perturbation theory. In Sect. 5 we discuss also the 
corresponding constraints on properties of the vortices. 

2. Perturbative expansions 

2.1. Divergences of perturbative expansions 

Standard analysis of divergences of perturbative expansions ^ introduces 
hierarchy of three scales. 



''This subsection is a very sketchy review. For details, see, e.g., 
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^QCD < « ^Iv . 

where Aqqd and Ajjv are the standard infrared and ultraviolet scales while 
the intermediate scale Q is specific for a particular process considered. Then 
a generic perturbative expansion for an observable ( O ) looks as: 

oo 

(O) = {parton model) ■ {1 + ^a„<((3^)) , (4) 

n=l 

where for simplicity of presentation we normalized the anomalous dimension 
of the operator O to zero. Note also that as{Q^) ^ 1. 

In fact, expansions (4) are only formal since the coefficients a„ grow 
factorially at large n: 

\an\ - c1-nl , (5) 

where Cj are constants. Moreover, there are a few sources of the growth 
(5) and, respectively, Cj can take on various values. First, the growth of 
the coefficients a„ may reflect the growth of the number of perturbative 
graphs in high orders, which is combinatorial in nature. The value of q 
is determined then by the classical action of an instanton solution. On 
the other hand, the divergence (5) can be triggered by a single graph of 
n-th order. This is a renormalon-type graph which contains n insertions of 
vacuum polarization into a gauge boson propagator. 

The factorial growth of o„ implies that the expansion (4) is asymptotic 
at best. Which means, in turn, that (4) cannot approximate a physical 
quantity to accuracy better than 

A ~ exp(-l/c,a,(Q^)) ~ ^^^y"^'' , (g) 

where bo is the first coefficient in the /3-function, 

Q^^asiQ') = -boaliQ') - b^aliQ') - ... , 

and we accounted for the fact that as{Q^) is logarithmically small at large 

To compensate for these intrinsic uncertainties one modifies the original 
expansion (4) by adding power corrections with unknown coefficients: 

oo 

(O) = {parton model) ■(!+ ^a„Q:(02) + ^bkiAgcD/Q)'' ), (7) 

n=l k 
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where powers k are determined in terms of cj. entering Eq. (5) and the fac- 
torial divergences are removed from the modified perturbative coefficients 
a„, for further discussion see, e.g., 

2.2. Borel summation 

So far, we discussed only behavior of absolute values of the coefficients a„. 
Now, if there is sign oscillation, 

a„~(-l)"n!cr , (8) 

then the sum is Borel summablc while if a„ ^ (+l)"n!c" there is no way 
at all to define the sum. In case of an asymptotically free theory the sign 
oscillation (8) is characteristic for the ultraviolet renormalons. For the first 
ultraviolet renormalon, 

(a„)„„ ~ n\bZ{-lT , (9) 

where bo is the first coefficient in the /3-function. Usually one does not 

reserve for any uncertainty in case of a Borel summable scries. 

The criterion of Borel summability might look too formal. Let us men- 
tion, therefore, that there exist also intuitive reasons to believe that there is 
no intrinsic uncertainty of perturbative scries due to the ultraviolet renor- 
malons. Indeed, one can readily check that the physical meaning of the 
power terms is that they correspond to contributions of non-typical vir- 
tual momenta k^^^^ which are either small or large compared to Q"^. The 
ultraviolet renormalons are associated with 

while infrared renormalons correspond to 

Furthermore, it is rather obvious that the contribution of large momenta 
in asymptotically free theories is calculable exactly. Then, there could be 
no intrinsic uncertainty due to the ultraviolet renormalon. 

On a more technical level, the argument runs as follows If one 
changes normalization of the coupling used in the expansion (4) from g^{Q^) 
to g^{p?) then the uncertainty (6) associated with the ultraviolet renor- 
malon changes to: 
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Thus using a normalization point fi^ 3> one can get rid of the corre- 
sponding power-hke uncertainty. There is no similar trick in case of the 
infrared renormalons, signaling that the corresponding power terms is a 
genuine uncertainty of the perturbative expansion. 

2.3. Non-perturbative match to infrared renormalon 

In phenomenological applications, we will concentrate mostly on two par- 
ticular examples, that is, gluon condensate (q;s(G'J^j^)^), where G"^ is the 
non-Abelian field-strength tensor, and the heavy quark potential V{r) at 
short distances, r — > 0. Keeping only the leading power corrections we have 

13. 

{{G%r) « {N^-l)-a-\l + ^ar.a:{a') + c^AQCD-a)'), (11) 
where Nc is the number of colors and a is the lattice spacing, and I4,i5,i6. 

limF(r) « + J2^^a^{r) + CsA^cd^^ ) . (12) 

A salient feature of the predictions (11) and (12) is the absence of quadratic 
corrections of order {Aqcd • a)^, {Aqcd ■ r)^- 

It is worth emphasizing that although the power corrections are so to say 
detected through pure perturbative graphs the actual vacuum fluctuations 
which dominate the power corrections in (11), (12) can well be genuinely 
non-perturbative. In particular, the leading power corrections both in (11) 
and (12) correspond to the so called infrared renormalon, a perturbative 
graph with many iterations of a loop insertion into the gluon propagator. 
However, the same dependence on Aqcd is provided by instantons: 

{as{G;,)%nst « {const)Al^cD ■ (13) 

Moreover, there are well developed and, in many respect, successful models 
of vacuum which assume the instanton dominance, see and references 
therein. 

It is important to realize that the instanton dominance in no way con- 
tradicts (11), (12). Indeed, analysis of the n-dependence of the coefficients 
a„ fixes the position of a singularity in the Borel plane and predicts in 
this way that the leading power correction is proportional to Aqqj^. This 
prediction is confirmed by evaluating the instanton contribution. Instan- 
tons might enhance the correction numerically and make the whole analysis 
more tractable. 
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2.4. Status of the ultraviolet renormalon 

Turning back to the ultraviolet renormalon, - which is most relevant to the 
present paper, - one might feel that the summary presented above is in 
fact self-contradictory. Indeed, if we start with expansion in UsiQ"^), see 
(4) then wc would reserve for a (Ag^^^^/Q^) correction corresponding to the 
ultraviolet renormalon. Which would be the leading power correction. Even 
if we apply summation a la Borel, a {^qcd/Q^) contribution would arise 
as a result of the summation of the divergent tail of the perturbative series. 
On the other hand, if we start with expanding in a^^^'^), fi^ 3> Q^, there is 
no intrinsic uncertainty to the perturbative series which would correspond 
to the ultraviolet renormalon. 

The contradiction is superficial, however. Namely, if one uses expansion 
in as{iJ^), then the {^qcd/Q^) terms are to emerge as a result of explicit 
summation of high orders of perturbation theory. Indeed, in terms of the 
original expansion (4) the ultraviolet renormalon corresponds to terms of 
order 

In terms of the as{n^) expansion one needs to keep even higher orders. 

It is worth emphasizing that an explicit calculation of the ultraviolet 
renormalon in gluodynamics is not straightforward at all The point is 
that the ultraviolet-renormalon divergence (9) is in fact related not only 
to the simplest chain graph but to a whole class of graphs, and no explicit 
calculation is possible. 

To summarize, one expects that the quadratic corrections are present 
but they are hidden either in the postulated procedure of summation a la 
Borel or in explicit summation of higher orders. It is worth emphasizing 
that numerical enhancement of the [Kq^^ jj / Q"^) terms is not ruled out by 
this analysis. However, such an enhancement is not suggested either. 

3. Status of corrections 

3.1. Power corrections from the infrared 

There exists huge literature on power corrections, for review see, e.g., 
Here, we will mention only a few points relevant to our discussion. 
The most standard sum rules apply to the two-point functions 

n,(Q2) = i /'exp(igar)(0|T{i(a;),i(0)}|0), (g^ = -Q'') , (14) 
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where j(x) is a current constructed on the quark and gluon fields and for 
simphcity we omitted the Lorenz indices. Moreover, one usually considers 
the Borel transform of (14), 

n,(Af2) = [Q^y{n-iy.](-d/dQ'rU,{Q') , 

where n — > oo, — > oo with Q'^/n = fixed. 

Somewhat symbolically, the sum rules then read: 

n,(M2) « (parton model)(l + + ^ + O {(HMy A^^ d))-'. 

(15) 

where the constants aj,Cj depend on the channel, i.e. on the current j. 
Terms of order 1/ In and are the first perturbative correction and 
contribution of the gluon condensate, respectively. Note that we kept only 
first order correction in a(M^). We could have included, say, one more per- 
turbative term, but the sum-rule approach is still based on the assumption 
that the effect of the power corrections is most important numerically. 

The physical meaning of the power corrections in (15) is that they pa- 
rameterize contribution of large distances of order Ag^^ in terms of local 
operators, like (GJ^^)^. The technical means behind is the operator product 
expansion (OPE) and there the procedure is well defined and unambiguous 

20 

While there is no doubt that (15) retains the leading power corrections 
if is large enough, it was noticed quite early "^^ that in some channels 
there could exist terms which start as sub-leading power corrections at large 
but numerically dominate at moderate ~ few GeV"^. Moreover, 
realistically such corrections are associated with so called direct instantons 

for a very recent and impressive example see 



3.2. The case for large Q ^ corrections 

While sum rules (15) claim many successes, there has also been accumulat- 
ing evidence suggesting introduction of Q~'^ corrections which apparently 
go beyond the standard picture summarized in the preceding subsection. 

To begin with, the instanton-dominated model of vacuum does not re- 
produce the confining potential for heavy external quarks On the other 
hand, lattice studies indicate that the so called Cornell potential for heavy 
quarks, 

^QqW « + cT-r , (16) 
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describes the data at all the distances r, for a review sec At large 
distances, the linear potential provides confinement. For our purposes here 
it is crucial, however, that the fit (16 suggests also that if we start with 
short distances the leading power correction to the Coulornb-likc intc^raction 
at 'moderate' distances r is of order a ■ r^. This observation is one of the 
motivations to modify the r.h.s. of Eq. (12) by adding a {Aqcd ■ term 
(or 5V{r)^Klcj,T) 

The gluon condensate, also, contains large quadratic corrections if the 
perturbative series is truncated at first terms. We will discuss this example 
in detail in the next subsection. 

Another dramatic example of large quadratic corrections is found in 
Ref. and concerns instanton density as function of the instanton size p. 
Namely, according to the standard operator product expansion, underlying 
also the sum rules (15), one has 

dn{p) « dn,ert{p){l + ^{0\{G;^m) , (17) 

where g{p) is the gauge coupling and the salient feature of (17) is, again, 
absence of a quadratic correction. However, it is demonstrated in that 
the data on the instanton density, to the contrary, unequivocally require 
a quadratic term. The data on the instanton density are described by an 
efii'ective action: 

g^2 

^^ff ~ ;^2(^ + '^P'^P'^ ' Cp ^ 2w , (18) 

where the first term is the standard action in the quasiclassical approxima- 
tion and a is the string tension. The expression (18) describes the variation 
of the instanton density dn{p) by about four orders of magnitude. 

There exist further examples when quadratic corrections improve phe- 
nomenological fits, see, e.g., However, the corrections are associ- 
ated with short distances and there is no substitution for the OPE which 
would allow to relate quadratic corrections to various quantities in a model- 
independent way. One of the ideas put forward in different contents, see, 
e.g., is a universal change of the running coupling by a Q~'^ term. In 
particular, one argues sometimes that analyticity requires removal of 
the Landau pole from the running coupling. The use of the modified, or 
'analytical' coupling, 
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then introduces a l/Q^ correction at large Q^. An apparent weak point 
is that imposing analyticity on higher orders in introduces further 
corrections which are not suppressed 

One can try to generaHzc (19) and think in terms of a universal coupling 
which includes a term. Let us note that the quadratic corrections in 
(16) and (18) differ numerically by a factor of 1.5 and agree in sign if the 
correction is ascribed to a universal coupling. This can be considered as 
a success of the model. Anyhow, till now there was no more convincing 
theoretical estimates of the quadratic term in (18). 

Quantitatively, the model with a short-distance gluon mass turns to 
be most successful. According to the model, the sum rules (15) are modified 
as 

n,(M^) . (parton model) (l + .^J/^^^^ + + +- ) ' ^''^ 

where the coefficients bj are proportional to the gluon mass squared and 
are calculable for any channel. The model cleared all the hurdles known. 
For example, one finds 

b„ « 4&P , (21) 

and this resolves a long standing puzzle of the analysis of the sum rules 
in the 7r-meson channel. 

To summarize, it seems reasonable, - as far as phenomenology of the 
confinement-related effects is concerned, - to replace the standard free gluon 
propagator by a propagator which reproduces the whole of the potential 
(16) already in the zero-order approximation of perturbation theory. On the 
theoretical side, such a program is a refinement of the variational approach 
•^^ and is outlined in Practically, full higher order calculations within 
such a scheme are known to be very cumbersome As far as the power 
corrections are concerned, this approach amounts to introducing a tachyonic 
gluon mass at short distances already in the lowest order of perturbation 
theory. Which is easy to implement and successful phenomenologically 

3.3. Perturbative - non-perturbative 'ambiguity' 

Note that all the phenomenological successes mentioned in the previous 
section refer to the fits with only first order perturbative contribution re- 
tained, along with a term. The successes claimed refer to moderate 
where the power correction becomes sizable. However, one could argue 
that if we start with large then to extract a small power-like corrections 
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one needs to subtract as many orders of perturbation theory as possible. 
Then, the question is how this would affect the fits to the correction 
at small distances. 

Numerically, this problem was studied in greatest detail in case of the 
gluon condensate on the lattice 35.36 j^^ terms of the lattice formulation 
the gluon condensate is nothing else but the average plaquette action . 
The result of the calculations can be summarized in the following way. 
Represent the plaquette action (P) as: 

(P) « P,^,, + h^'a'KlcD +c''a^^'QCD , (22) 

where the average plaquette action (P) is measurable directly on the lat- 
tice and is known to high accuracy, Pj^^^ is the perturbative contribution 
calculated up to order N: 

n=N 
n=l 

and, finally coefficients are fitting parameters whose value depends 

on the number of loops N. Moreover, the form of the fitting function 
(22) is rather suggested by the data than imposed because of theoretical 
considerations. 

The conclusion is that up to ten loops, TV = 10 it is the quadratic 
correction which is seen on the plots while are consistent with zero. 
However, the value of decreases monotonically with growing N '^^ and 
somewhere after 10 loops becomes consistent with zero while the A.qqj-, 
term finally shows up Moreover, the emerging value of the Aq^j^, term 
is consistent with current phenomenological estimates from the sum rules 

36 

Turn now to another example of a large Q~'^ correction, that is the 

heavy quark potential, see (16). The first attempt to analyze the effect of 
higher perturbative corrections on the linear term at short distances was 
undertaken in Ref. Defining the potential as 

^gg(«) - -C.l^ , (24) 

and using first terms in the perturbative expansion of ay which are known 
explicitly, 

aviq) = aMsiq){^ + ai^MsiQ) + (^2a\jg{q) ) , (25) 

one finds that the linear piece in the potential at short distances is 5 times 
larger than it would follow from (16). In other words, the result depends 
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crucially on the subtraction procedure. To proceed further, one is to invoke 
a model. 

In particular, it was suggested to saturate higher orders by the leading 
infrared renormalon evaluated in the large-6o approximation: 

N 

av{q) ~ aM5(g)(l + aia^s + O'^a^Ms + Yl°'n''oiMs) (26) 

where N ~ 3/26oQ!ms- The potential observed on the lattice is reproduced 
then by (26) including 5V{r) ^ ar at short distances. Moreover, the re- 
sult is stable against reasonable variations in N. Note that no explicit 
terms are to be introduced within this procedure since high orders of 
perturbation theory are presumably accounted for explicitly. 



3.4. Complementary ways of describing the Q ^ terms 

Thus, the lesson from calculations of the perturbative gluon condensate is 
that one can either approximate total matrix elements by a low-order per- 
turbative contribution plus a quadratic correction, or by high-order pertur- 
bative contributions plus a quartic correction. 

Our central point, which motivated us to review the evidence in favor 
and against non-standard corrections, is that there is no contradiction 
between the two approaches. Moreover, existence of the dual descriptions of 
the terms is expected in fact on pure theoretical grounds, see discussion 
of the status of the ultraviolet renormalon in Sect. 2.4 . 

Theoretically, it is known that no ad hoc terms are allowed to be 
added to untruncated perturbative series However, then one is to be 
prepared to calculate high orders indeed. In the only case when such a 
calculation turns possible (that is, the gluon condensate ^5,36^ theoretical 
expectations are fully confirmed. Namely, lowest orders plus a quadratic 
correction give reasonable fits 35,36 ^j^g other hand, if many loops are 
accounted for, then there is no need for an ad hoc power correction any 
longer What remains unanswered at this point is why the terms, - 
described in one or the other way,- are important phenomenologically. As 
we will argue later, observation of the branes seems to provide us with a 
key to answer this question. 
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4. Branes and power corrections 

4.1. Quadratic correction to the gluon condensate 

Monopoles and vortices are detected, for a given configuration of the vac- 
uum fields, for the whole of the lattice. Thus, they are seen as a nonlocal 
structure. Moreover, they are manifestly non-perturbative. Indeed, the 
probability 9{plaq) for a particular plaquette to belong to a brane has been 
found to be proportional to: 

9{plaq) Rs {cons t)exp( -l/bog\a)) ^ (a ■ Aqcd)' • (27) 

On the other hand, the branes have an ultraviolet divergent tension which 
assumes a kind of locality. 

To make contact with the continuum theory it is useful to evaluate 
contribution of the branes into local or quasi-local matrix elements. The 
gluon condensate (11) turns to be the easiest case. Indeed, combining 
Eqs (3) and (1) one gets for the contribution of the vortices to the gluon 
condensate: 

< (G^.)' )vort ~ 0.3 GeV^ a-^ , (28) 

which matches the ultraviolet renormalon. The beauty of this result is that 
all the quantities considered are manifestly gauge invariant. 

4.2. Quadratic correction to the heavy quark potential 

In case of the heavy quark potential, one can also argue that the branes 

match the ultraviolet renormalon, or the correction. Indeed, it is well 
known that both monopoles and vortices generate the linear potential for 
heavy quarks. In particular, in case of the monopoles *°: 

(^QC3M)™o„ « ''rnon-T , (29) 

where a is close numerically to the string tension in the full potential (16). 
Moreover, the potential (29) is linear at all distances tested, beginning with 
r = a. The reservation is that the observation (29) refers to the Abelian 
projected potential. However, the approximation is known to be valid for 
the quarks in the fundamental representation ( for detailed discussion see, 
e.g., ^). 

4.3. Probing properties of the branes 

To produce a linear potential (29) at short distances, vacuum fluctuations 
are to satisfy highly non-trivial constraints. First, the size of fluctuations is 
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to be small i^'^i. Indeed consider again for simplicity the Abelian projec- 
tion. Then the potential can be represented as an integral over the overlap 
of electric fields associated with the external charges charges ±Qei- 

V{r) = E(r')-E(r + r')dV , (30) 

where r is the radius-vector connecting the charges. Imagine that there 
is a change in the field at distances of order R. In the applications, R ~ 
{Aqcd)~^- At distances R^ r the field of the charges is a dipole field: 

lEd.pl - Qr/R^ . 

Respectively, the change in the potential is of the order: 

SV{r) ~ l^d^r'/ir'f ~ ^ , (31) 

and we would conclude that the non-perturbative potential at small dis- 
tances is proportional to r^, see also Eq. (12). 

Therefore, the standard estimate might not work only in case that fluc- 
tuations responsible for the confinement have size comparable or less than 
the distance between the quarks, r. Monopoles and vortices with the action 
(1) do satisfy this constraint since the ultraviolet divergence in the action 
implies size of order a. 

To appreciate another constraint on the properties of the vortices, con- 
sider the Z2 projection, (for a review see ^). Then the area law for the Wil- 
son loop is derived from the assumption that the probabilities for plaqucttes 
in the plane of the Wilson loop to belong to the vortex are uncorrelated 
Indeed, in this case the Wilson loop on average is given by: 

{W{C)) = Ilpia,{{l-f) + .n-l)){WoiC)) , (32) 

where the product is taken over all plaquettes belonging to the minimal area 
A spanned on the contour C, / is the probability for a plaquette to belong 
to the vortex, and {Wo{C)) is the value of the loop with the constraint that 
no vortices pierce the minimal area. Then, 

{WiC)) « exp{-aA/a^){Wo{C)) , (33) 

where the string tension cr = — ln(l — 2/) and A is the minimal area in 
the lattice units. Note that / ^ (aAgczj)^; see Eq. (27). 

Now, we are discussing the case when the contour C has one of dimen- 
sions very small, A T • r where T is the time extension and r ^ 0. Eq. 
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(32) applies only if the vortex remains random even at distances compara- 
ble to the lattice spacing a. Which implies, in turn, that the entropy factor 
for the vortex grows exponentially with its area: 

{Entropy)yort ~ exp( cA/a^) . (34) 

This is indeed true for the branes observed on the lattice, see ^. 

So far, we have not discussed corrections to the instanton density, see 
(18) since no quantitative framework is known in this case. Let us only 
note that the physics behind the standard prediction (17) is similar to the 
case of the heavy quark potential. Indeed, in case of instantons we deal 
with color dipoles in d=4. Again, the OPE prediction (17) follows from 
the assumption that the confining fields are soft and modify the original 
instanton field at distances of order {Aqcd)~^- The validity of (18) implies, 
therefore, that there exist small-size vacuum fluctuations. There are no 
other (non-perturbative) candidates but the branes. 



4.4. Branes and perturbative series 

The original ultraviolet renormalon sequence of the perturbative coeflicients 
is given by (9). The n\ factor arises then from integrals over virtual mo- 
menta of the type: 

y.g. ki • 

Such a contribution is related to a single graph with n insertions of the 
vacuum polarization into a gauge boson propagator. However, this is not 
a single source of contributions of order (9) Namely, one can reserve, 
say, for one loop insertion less because then one gains in terms of number 
of graphs. In this way there arise contributions of order 

an ~ {n — 1)\ ■ I , I ^ n , 

where the first factor is due the mechanism (35) and the factor Z ~ n is 

combinatorial. Thus, one gradually switches to evaluating both the number 
of graphs and powers of the logs. This two-parameter problem becomes 
practically intractable after a few steps. 

It is worth emphasizing that in case of the power corrections to the 
gluon condensate there is no mechanism (35) at all. Indeed the expansion 
is in terms of g'^{Afjy) and there are no virtual momenta larger than the 
normalization point. 
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Thus, the only mechanism left for having a large quadratic correction 
is a large number of graphs. Moreover, large power-like terms cannot be 
associated now with sign oscillation. Indeed, sign oscillation plus expansion 
in (7^(A^y) would result in a negligible contribution. Thus, we come to 
conclusion that the series associated with the quadratic correction should 
not oscillate in sign. Only then we can expect that the series corresponds 
to a sizable quadratic correction. 

Let us check the expectations against the only known example of an 
explicit calculation of a 'long' perturbative series, that is, turn to the case 
of the gluon condensate in the lattice SU{3) gluodynamics. According to 



36. 



~ u(l - ^) , (36) 



Pn-l V n+s 

where the coefficients Pn are introduced in (23) the numerical values of the 
parameters are: 

u = 0.961(9), q = 0.99(7) , s = 0.44(10) . 

The quadratic correction is affected by 10. 

We see, indeed, that the perturbative series has no sign oscillation. 

So far, two regular mechanisms for generating same-sign perturba- 
tive expansions were discussed Namely, infrared rcnormalons and 
perturbative-vacuum instability due to instantons. In both cases the ra- 
tio r„ grows at large n linearly with n. In case of the leading infrared 
renormalon: 

r-r - ^ • (37) 

Remarkably enough, in the crucial region of n ~ 10 this ratio is still smaller 
than (36) and the infrared renormalon catches up with (36) only around 
n ~ 25 In case of instanton-related divergence the ratio r„ is even 
smaller than in case (37): 

rr* « ^ . (38) 
" 447r ^ ' 

Prom the theoretical point of view, both the instantons and infrared 

renormalons are irrelevant to the corrections. Thus, it is gratifying 
that the perturbative series (36) looks indeed different. The series (36) 
is convergent for \g'^\ < u and seems to exhibit a novel mechanism of 
generating same-sign perturbative series. 
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4.5. Crossover 

The non-analyticity in g"^ indicated by (36) might be a reflection of the 
crossover transition Indeed, the irregular behavior of the specific heat 
associated with the crossover is known since long Note also that at 
the crossover the brancs, with the properties known in the wcak-couphng 
region, disappear. Indeed, in the strong-coupling region, that is beyond 
the crossover there is no scale Aqcd at all and the basic properties of the 
branes, like (27) make no sense. 

Thus, there arises the following tentative picture of the mechanism be- 
hind the corrections. Because of the crossover transition, there appear 
corrections to the standard running of the coupling at mass scale of about 
GeV. This change in the running is manifested phenomcnologically in var- 
ious ways. Thus, there is evidence for a 'freezing' of the effective coupling, 
for review see, e.g., In particular, the following effective coupling: 

bolmod ^ ^0 ^mod ' 



where 



^■mod ^ m - 



^ ^^QCD 



fits well various pieces of data ^"^ . 

On the non-perturbative side, the branes with properties (1), (3) emerge 
around the crossover (if one moves towards weak coupling). The simplest 
and qualitative fit is therefore lowest order in perturbation theory plus ef- 
fect of the branes. This fit is justified in the region close to the crossover - 
and mostly the lattice measurements are performed for (3 not far from the 
crossover. On the other hand, the same effect is to be described by per- 
turbation theory since the measurements are already in the weak-coupling 
domain and there is no singularity which would block the pcrturbative ex- 
pansion from being valid in the region close to the crossover (approached 
from the weak-coupling side). The price is that, the closer is the crossover, 
the more terms in the pcrturbative expansion are to be kept. Evaluating 
so many pcrturbative terms is practically impossible except for the case of 
the gluon condensate 

The pcrturbative series corresponding to the Q"^ correction exhibits 
sign coherence, similar to the case of a classical solution. The series is not 
divergent in large orders, however. That is, branes are no classical solution. 
Numerically, the pcrturbative expansion coefficients in the crucial region of 
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n are even larger than for known divergent series (infrared renormalon, 
instantons). This is correlated with the fact that the branes apparently 
represent strongly aligned vacuum fluctuations. 

5. Constraints on the branes 

5.1. Consistency of the branes with asymptotic freedom 

In the preceding section we could convince ourselves that to match the 
ultraviolet renormalon the branes are to have highly non-trivial proper- 
ties like (1), (27), (34). All these properties, crucial to fit the ultraviolet 
renormalon, reveal a point-like facet of the branes. Now we arc revers- 
ing the question and ask whether this point-likeness is consistent with the 
asymptotic freedom. In particular, the action for the lattice monopoles is 
proportional to a~^, the same as for point-like particles. At first sight, it is 
Tinavoidable that accepting (1) in the limit a ^ is equivalent to introduc- 
ing new particles at short distances. Appearance of such particles would 
be inconsistent with the asymptotic freedom. 

To have a closer look at the problem it is useful to translate the data on 
the monopole trajectories into a conventional field theoretic language. To 
describe monopoles one then introduces a magnetically charged field 0m- 
Moreover, since the monopoles are observed as trajectories, it is natural 
to use the polymer representation of field theory, see, in particular, 
Proceeding in this way one can derive (see and references therein): 

(0| |(/.m|' |0) = ^{pperc + Pf^n) , (40) 

where density of the percolating cluster, Pperc is defined in (2) while pfin 
is related in a similar way to the length of the finite monopole clusters. 

The total density of the monopole clusters can be measured directly 
It is crucial that the total monopole density diverges as at small a: 

Pperc + Pfin ~ {const)hl^cD + {cOUSt)—^ . (41) 

Finally, we get: 

(0| \cbM\^ |0) « [const )AIcd/» « 0.8 (/m)-^. (42) 

Thus, (0| |0Mp |0) contains no ultraviolet divergence and is, therefore, 
perfectly consistent with the asymptotic freedom which does not allow to 
add new particles. 

To appreciate the geometrical meaning of the observation (41) turn 
to the simplest case of uncorrelated percolation. One introduces then a 
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probability p for a link to be "open", that is to belong to a monopole 
trajectory in our notations. Then at a critical value pcr there arises an 
infinite, or percolating cluster. The density of this percolating cluster is 
vanishing at the point of the phase transition to the percolation. In the 
supercritical phase where p > Pcr this density 

Pperc ~ {p-PcrT , (43) 

where the critical exponent a > 0. 

While the density (43) is vanishing at p = Pcr, no non-analyticity can 
happen in the total density. At the point of the phase transition, it is given 
simply by: 

Ptot = ^ , (44) 

where d is the number of dimensions. 

Substituting (44) into (40) we recover in case of d = 4 the standard 
quadratic divergence, 

(0| \ct>M\' |0)~a-2 

which is so familiar from the case of Higgs particles. However, the data 

(41) indicate only behavior of the total monopole density. Which geo- 
metrically means that the monopoles spread over a = 2 subspace of the 
total d = 2 space. (Of course, the d = 2 subspace is not a plane but rather 
a surface percolating itself through the d = 4 space). 

Thus, we come to an amusing conclusion that it is the existence of the 
branes which eliminates a potential quadratic divergence in (40). Note 
that the fact that the monopoles are associated with the vortices (whose 
total area scales) was observed first in Ref. Later the phenomenon was 
confirmed for various values of a In view of the ultraviolet divergence 
in the monopole action, see (1), this association of the monopoles with 
vortices becomes absolutely crucial for the consistency with the asymptotic 
freedom. 

Our final comment concerning (42) is that the vacuum expectation value 

(42) is perfectly gauge invariant. Gauge invariant condensates of dimension 
two were widely discussed recently, see, e.g. The beauty of the relation 
(42) is that it does not contain ultraviolet divergences which plague local 
condensates formulated in terms of the gluon fields. Thus, the vacuum 
expectation value (42) could be related directly to the non-perturbative 
part of the gluonic dimension two condensates. No explicit relation of 
this kind is known, however. 
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5.2. Matching of the branes with ultraviolet renormalon 

The lattice data on the branes, see, in particular (1), (3) are compelling 
but pure numerical. Since structure of the branes in terms of the original 
gluon fields is not known, it is very difficult, for example, to explain (27) 
theoretically. Matching of the branes with the ultraviolet renormalon brings 
constraints on the branes which turn to be quite restrictive, fn particular, 
to match the ultraviolet renormalon in case of the gluon condensate the 
branes are to obviously satisfy the following constraint: 

^ (thickness) /, , ■ ^ 4\ .2 2 /.r.^ 

f'piaq ■ ■ ((plaquctte action) ■ a ) ~ Ag^j^a , (45) 

where 9piaq is the probability for a given plaquette to belong to the vortex 
(see (27)) and the thickness is understood in terms of the distribution of 
the excess of the SU{2) action. The presently available data indicate ^ that 
0piaq ~ i^QCD ' d)^ , the plaquettc action is of order and the thickness 
is equal to the resolution, that is a. (All the quantities are on the average.) 
Such a regime is fully consistent with (45) and can, therefore, persist in the 
limit a ^ as well. 

There is another and subtler point. The non-Abelian action, albeit 
ultraviolet divergent, is known to be much smaller than the projected one. 
In case of monopoles this was emphasized, in particular, in 

Qmon ^ , ^ Qmon _ ^ L / 

a a 

where 5*^';,°" is the action in the Abelian projection, i.e. the same as for the 
Dirac monopole. Similar inequality holds in case of the vortices. 

Now, wo can dcirive (46). Indeed, the })ranes are dual to high orders 
of perturbation theory. Which moans that the action associated with 
the branes is is not allowed to be parametrically more singular than the 
perturbative fluctuations. This rules out the monopole action of order 
Smon ^ ■ L/a, see Eq. (46). In other words, branes cannot be removed 
from the vacuum without affecting perturbative fields as well. 

5.3. Casimir scaling 

The lattice data on the quark potential (for review see, e.g., '^^) reveal a 
strikingly simple picture. Lack of a same simple theoretical explanation 
looks as a puzzle. In fact there exist rather a few puzzles to be explained: 
a) at relatively short distances the potential is well approximated by a 
Coulomb like piece plus a linear potential; 
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b) the linear piece continues, with the same slope to large distances; 

c) Casimir scaling. 

The phenomenon of the Casimir scaling is that at intermediate dis- 
tances static potentials for quarks belonging to various representations D 
of the color group are proportional to each other: 

Vn{r) « p^VV(r) , (47) 

where Cd is the eigenvalue of the quadratic Casimir operator 

Cd = TrT^T^ 

of the representation, Cp is the value oi Cn for the fundamental represen- 
tation and Vpir) is the static quark potential in case of heavy quarks in 
the fundamental representation. 

Explanation of the observed Casimir scaling is a challenge to theory, for 
review and models see 50,5i^ obvious either that observation of the 

branes brings a transparent solution to the problem. We still feel that the 
consideration of the branes does introduce new elements into understanding 
of the Casimir scaling. 

As far as the observation b) above is concerned, a new point is that the 
branes look the same random at all scales tested. Namely, the properties 
(1), (3), (34) reiterate themselves for all the lattice spacings. Since there 
is no scale involved, except for Aqcd, see Eq. (27) the same slope a ~ 
Aqcd governs the brane-induced potential at all the distances. An explicit 
realization of this idea is provided by Eq. (32). 

Thus, the branes allow to shift the emphasis in explaining the Casimir 
scaling from large to short distances. Concerning short distances, the linear 
potential is calculable, in principle, perturbatively. In the two-loop approx- 
imation (that is, three terms in the expansion of V{r) in as) the Casimir 
scaling is checked by explicit calculations 

ytwoioop^^^ = g£v^;'"°-'°°P(r) . (48) 

However, the linear term at short distances is sensitive to even higher orders 

(for a discussion see Sect 3.3) and one has to turn to models. In particular, 
the model with a tachyonic gluon mass immediately gives: 

SVnir) = ^\ml\r , (49) 



reproducing the Casimir scaling. 
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More generally, the branes are to be dual to the ultraviolet renormalons. 
In turn, the factorization properties of the ultraviolet renormalon, see 
and references therein, result in the Casimir scaling. However, within such 
a framework it is difficult to expect that the Casimir scaling would hold to 
high accuracy. 

To summarize, consideration of branes allows to reduce the problem of 
the Casimir scaling to the problem of evaluating the linear potential at short 
distances. The latter problem is perturbative in nature and the first two 
loop corrections are known to exhibit the Casimir scaling. Higher orders 
can be estimated only within models. 

6. Conclusions 

We have demonstrated that the monopoles and vortices are responsible for 
the ultraviolet-renormalon type corrections at short distances. In turn, the 
Q~'^ corrections are dual to high orders of perturbation theory. Phenomeno- 
logically, such corrections are known to be welcome. Combination of the 
ultraviolet and infrared factors typical for the ultraviolet renormalon ap- 
pears to be a reflection of the fine tuning for the monopoles and vortices as 
non-local objects. Vortices and monopoles appear, in this context, as non- 
perturbativc counterpart of the ultraviolet renormalon. Highly non-trivial 
constraints on the properties of the branes implied by this identification are 
satisfied by the lattice data. 

The overall conclusion is that the fine tuning observed on presently 
available lattices, see in particular (1), (3), could be the true asymptotic in 
the limit a — *■ 0. 

In a broader context, fluctuations appearing as topological and sup- 
pressed in one formulation of a theory can become fundamental entities 
in a dual formulation of the same theory. If it is true in case of SU{2) 
Yang-Mills theory, then the hint is that the branes could appear in the 
dual formulation. Where by branes we understand (see above) d = 2 sur- 
faces populated by the monopoles (tachyonic mode) and living on a d = 4 
Euclidean space. Generically, this observation agrees with recent and well 
known proposals on theories dual to YM theories. To the best of our 
knowledge, no direct comparison of the two approaches is possible, however. 
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